Diffusion of hydrogen in Pd assisted by inelastic ballistic hot electrons 
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Sykes et al. [Proc. Natl. Acad. Sci. 102, 17907 (2005)] have reported how electrons injected from 
a scanning tunneling microscope modify the diffusion rates of H buried beneath Pd(lll). A key 
point in that experiment is the symmetry between positive and negative voltages for H extraction, 
which is difficult to explain in view of the large asymmetry in Pd between the electron and hole 
densities of states. Combining concepts from the theory of ballistic electron microscopy and electron- 
phonon scattering we show that H diffusion is driven by the s-band electrons only, which explains 
the observed symmetry. 

PACS numbers: 63.20.kd, 68.37.Ef, 71.20.Be 
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While the scanning tunneling microscope (STM) seems 
a natural tool to manipulate all sort of atoms and 
molecules adsorbed on surfaces [1-4 , it is rather surpris- 
ing to learn that tunneling currents can also be used to 
change diffusion rates of atoms deeply absorbed in met- 
als. In an elegant experiment performed on H buried be- 
low the Pd(lll) surface, Sykes et al. have unequivocally 
proved that this effect does exist [5 . Indeed, the STM 
has been used to investigate buried interfaces [6 or even 
individual molecules [7 , but the pioneering experiment of 
Sykes et al. confirms the STM as a powerful nanometric 
technique to control diffusion and reaction processes at 
the atomic level not only on surfaces but also inside met- 
als. The controlled insertion and extraction of H in Pd is 
an added interest of this particular experiment because 
of its possible applications in H storage [SHTO] and in 
heterogenous catalysis. In this respect, Pd is a recurring 
catalyst in hydrocarbon synthesis reactions, particularly 
in -C=C- bond breaking, and the reason behind such ex- 
ceptional catalycity is precisely linked to the existence of 
weakly bound H at subsurface sites [IJj. 

Still, the underlying mechanisms behind the selective 
H population of subsurface sites reported in [5 are not 
well understood. In particular, the observed bias volt- 
age dependence is puzzling. Working at low tempera- 
ture to prevent thermal diffusion (T = 4 K), the STM 
is used as a local electron gun to inject small currents 
{It ^ 1 — 150 p A) in the vicinity of the Fermi energy 
{V ^ ±1 V). Inelastic excitation with the injected ballis- 
tic hot electrons promotes the diffusion of H from bulk Pd 
towards the surface even under low concentration condi- 
tions, thereby implying the contribution of large volumes 
of the metal where electric fields are efficiently screened 
out and cannot play a role. The distortions in the Pd-Pd 
distances caused by subsurface H accumulation are then 
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FIG. 1. (a) Calculated density of states for bulk Pd. (b) 
Sketch for injection of s-like (isotropic propagation), and d- 
like (narrowly focused) electrons on Pd(lll). (c) Energy iso- 
surface {Ef + 0.83 eV) related to the s band in Pd. 



imaged with STM as protrusions or bright stripes. The 
brightness and width of the stripes increase with It and V 
magnitude. The intriguing finding is that the stripes are 
comparable in brightness for opposite V signs, positive 
biases creating only slightly wider stripes. This polar- 
ity effect must have a non trivial explanation, since for 
V > and V < charge carriers are electrons and holes 
respectively, and the density of states above and below 
the Fermi level Ep is very different in Pd [Fig. fl] (a)]. 

In the present letter, we examine the vibrational exci- 
tation of H in Pd caused by inelastic interaction with bal- 
listic hot electrons using first order perturbation theory 
for the electron-phonon (e-ph) coupling. Hopping rates 
for H migration have been quantified as a function of the 
experimental V and It. Making use of the theory devel- 
oped to describe propagation in ballistic electron energy 
microscopy (BEEM) [12 , we find that the observed sym- 
metry between positive and negative V follows when the 
different roles of s and d-band electrons has been recog- 
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FIG. 2. (Color online) Profile of the computed energy bar- 
rier (blue thick line and circles) for diffusion from the octahe- 
dral (O) to the tetrahedral site (T). The approximated har- 
monic potential in O and T (black thin lines) and energy levels 
(dashed lines) are given for comparison. When H is excited 
up to a vibrational level n > 2, transmission by quantum tun- 
neling becomes efficient and dominates in the intensity range 
of 1-150 pA. Inset: sketch of O and T sites (white) in the Pd 
fee lattice (dark). Arrows mark the diffusion path from O to 
T. 



nized in the H vibrational excitation. It has been noticed 
that d-like electrons in hulk Pd propagate along narrow 
cones corresponding to a few selected k points in the first 
Brillouin zone (IBZ) [13 . Such a focusing effect happens 
typically when, with increasing energy, bands approach- 
ing the borders of the IBZ change their curvature from 
convex to concave in order to touch the boundaries at 
right angles. While in a typical BEEM experiment fo- 
cused beams are used to improve the resolution well be- 
low the surface, these electrons span a reduced substrate 
volume and have a small probability to interact with in- 
terstitial H in the low dilution regime. On the contrary, 
electrons propagating through s channels span large vol- 
umes and are relevant for the present experiment because 
they are the only ones that can interact with many in- 
terstitial H [Fig. [1] (b) and (c)]. Interestingly our result 
shows that the inelastic excitation of H in Pd yields infor- 
mation about a component of the current that is hidden 
in a typical BEEM experiment, since it distributes over 
large spatial regions and does not contribute to nano- 
metric resolution. Furthermore, we find that at T = 4K 
the efficiency of the diffusion mechanism can only be un- 
derstood by incorporating quantum tunneling of H near 
the top of the barrier. Similar phenomena have been 
in fact observed for heavier elements, e.g., deuterium in 
[M]. For this particular case, a larger mass would reduce 
the tunneling rate, change the frequencies of vibrational 
modes, and modify the rates in general. However, the 
relevant physics shall remain the same. 



Ab-initio calculations have been performed by expand- 
ing wave functions on a plane-wave basis set using the 
Quantum Espresso package [15 . This formalism yields 
accurate total energies, equilibrium geometries, electron 
band structures and transition states, as well as phonon 
frequencies, eigenvectors (in first order density functional 
perturbation theory), and e-ph coupling matrix elements. 
These are the ingredients for the inelastic model de- 
scribed below. The actual experiment is mimicked by 
using a periodic supercell with PdieH stoichiometry. We 
have checked that such a supercell brings only small cor- 
relations between H atoms in consecutive cells. The the- 
ory level for the exchange and correlation functional is 
the generalized gradient approximation in the Perdew- 
Burke-Ernzerhof formulation [16 , and ion cores are de- 
scribed by ultrasoft pseudopotentials [17]. In the con- 
struction of the plane wave basis sets we have used 
Monkhorst-Pack /c-point meshes [18 and energy cut-offs 
for wave functions and charge densities of 28 and 180 Ry, 
respectively. The convergence threshold for total energy 
self-consistency is 10~^^Ry. The equilibrium lattice pa- 
rameter obtained for fee Pd is 3.98 A. Using a 14 x 14 x 14 
Monkhorst-Pack mesh and the PdieH cell, we determine 
that H is more stable at the octahedral site (O) than at 
the tetrahedral site (T) by 0.064 eV. The atomic posi- 
tions of the Pd atoms surrounding the H atom have been 
relaxed with tolerances of 10~^ Ry and 10~^ Ry/a.u. in 
the energies and forces, respectively. The energy barrier 
from O to the metastable T has been found by the nudge 
elastic band method [19 to be E^ = 0.178 eV (see Fig.[^. 
This is in good agreement with the 0.23 eV barrier found 
in electromigration experiments [20 . In the nudge elas- 
tic band calculation, nine images have been used and the 
final forces at the transition state are ^ 5 x 10~^ Ry/a.u. 

Phonons have been calculated at the F point using 
10~^^Ry as convergence criterion in the self-consistent 
loop. By restricting ourselves to q = modes, we are 
neglecting correlations between H atoms in consecutive 
unit cells, which is exact in the dilute limit. Only the 
H atom and its six nearest neighboring Pd atoms are 
allowed to move in the phonon calculations. The three 
nondegenerate vibrational mode eigenvectors that yield 
a significant displacement of the H atom have energies of 
17, 21 and 47meV. The constraint on Pd atoms affects 
these values by only ~ 5 meV. 

In order to calculate H diffusion from O to T sites we 
approximate the O well by a harmonic potential trun- 
cated at an energy equal to Ejj [2TH23] . Figure |2] shows 
that the highest frequency mode huj = 47meV, in good 
agreement with the barrier profile, is the most likely to be 
involved in this process. Therefore, the H atom ought to 
be vibrationally excited up to approximately the fourth 
harmonic level to overcome the barrier and diffuse. The 
excited level populations obey a master equation that, in 
the limit of weak inelasticity, yields the following transfer 



rate from the n th vibrational excited level 



Rn = nT 



ri 



(1) 



where T^ and F^ are the vibrational excitation and deex- 
citation rates between the ground and first excited states. 
Hence, Rn is obtained as the population of level n — 1 
times the probability to reach the nth vibrational excited 
state that is already above the top of the barrier. 

Hot electrons injected from the STM tip into the sub- 
strate, or in the inverse direction, have energies in excess 
according to the bias voltage, F, and can excite vibra- 
tional quanta {huj) on the interstitial H via e-ph cou- 
pling. In first order time-dependent perturbation theory 
the transition rates for electron {V > 0) or hole {V < 0) 
decay can be obtained as: 



rl'%;-) = ¥^E E 
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where k and g are the electron and phonon wavevectors, 
respectively. The t, | symbols account for n ^ n + 1 (vi- 
brational excitation) and n ^ n — 1 (vibrational deexcita- 
tion) transitions, respectively. The initial and final elec- 
tron band energies Cij lie in the range Ef < €ij < Ef-\-V 
for bias voltages F > 0, and Ef — \V\ < Cij < Ef for 
V < 0. The e-ph coupling matrix elements Mij{k^q;uj) 
are computed by Quantum Espresso [15 . Convergence 
on the sum over the IBZ in Eq. [2] has been checked by 
using a sufficiently fine /c-point grid [24 . A/'pd is the num- 
ber of Pd atoms in the unit cell, a normalization factor 
ensuring that the rates do not depend on the cell size. 

Now the crucial point is how to single out the injected 
electrons that are actually responsible for the H vibra- 
tional (de)excitation. With this purpose, we turn to the 
theory developed to describe propagation of hot electrons 
through a periodic lattice in BEEM [12]. In its simplest 
version, BEEM theory uses a semiclassical Greens func- 
tion due to Koster [25] to link the topography of the 
band-structure constant energy surface at Ef ± V with 
the spatial propagation of electrons in the substrate peri- 
odic lattice. For Pd, it predicts the formation of narrowly 
focused beams associated with the nearly flat d bands, 
while electrons injected in s bands propagate in a semi- 
spherical, 5- wave-like form (see Fig. fl]). For the experi- 
mental samples of diluted H in Pd we expect a behavior of 
the s and d bands similar to that found in bulk Pd, since 
the low H concentration has a small impact on the Pd 
band structure. This has been seen in the PdieH super- 
cell, that already represents well the low-density regime 
in the experiment [26]. Therefore, electrons injected in 
the strongly focussed d bands are not likely to encounter 
and, hence, to excite interstitial H. This mechanism can 
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FIG. 3. Bias voltage dependence of the vibrational excitation 
(rj) and deexcitation (FJ) rates as calculated from Eq. ([2l 
when including the propagation correction through w^.p (see 

ik 

text). Inset: FJ'^ if both the s and d bands contribute in the 
injected electrons. 



only be driven by the s bands propagating in all direc- 
tions in the crystal. To take into account this effect, 
we introduce in Eq. 2 weights, '^2% that have been ob- 
tained as the projected 5-orbital percentage in the whole 
\ilji{k)) Kohn-Sham state. The calculated FJ'^ are plot- 
ted in Fig. [Sl As seen in the inset, if both s and d states 

•—I >s I 

contributed fully, Fj'^ would map the Pd DOS and the 
large asymmetry around Ef would cause differences in 
the rates of at least 2 orders of magnitude between pos- 
itive and negative V. These differences are drastically 
reduced as soon as we realize that H excitation is mainly 
caused by nearly free (s-like) electrons, that have a sim- 
ilar dispersion above and below Ef- 

Rates in Eq.|2]are not yet the ones appearing in Eq.[T] 
The latter are the result of the net vibrational excitation 
caused by all the charge carriers in the incoming current. 
It. Therefore, we introduce the following linear model 
for F^'^: 



F^ 



It 

eTop 



F^ 
J- 1 5 



Fo 



It 

eVop 



^ 1 



(3) 



where Fq is the deexcitation rate present at zero bias [Fq 
is calculated from Eq. 2 simply by removing the w t 
factors], and p (see inset in Fig. [4|) is the number of 
available one-electron states per Pd atom. By inserting 
this factor we are assuming that all the incoming carri- 
ers have an initial energy e^ with equal probability 1/p, 
and also that the final state is always unoccupied. The 
latter is a fair assumption in the present low-intensity 
regime. Considering that s carriers distribute evenly 
throughout the crystal and that the typical time that 
a H atom spends in a given excited vibrational state is 
~ 1/Fo ~ 10~^^s, the factor ^ represents the fraction 



of the current that (de)excites the H atom. The weak 
inelastic hmit (F^ c^ Fq ^ F^) holds well for the ex- 
perimental intensities It = 1 — 150 pA. In particular, for 

V = ±0.7 V and It = 150 pA the effect of the injected hot 
electrons can be understood as if the interstitial H was 
effectively surrounded by a thermal bath at Teff ::^ 80K 
according to the expression for a Botzmann distribution 
of level populations A^^Teff = huj/\n{T^/T^) [22 . 

Quantum tunneling becomes important in the diffu- 
sion of nonmassive elements, like H, specially at low tem- 
peratures. We introduce it in the model via a semi- 
classical (WKB) transmission rate through the barrier, 
T{En) = e~'^^^^'^\ To a very good approximation, 
a triangular-shaped barrier can be used instead of the 
upper part of the pathway in Fig. |2j In this case, 

7(E,) = I [^{E, - En)]'^^d{Er,), where d{En) is the 
barrier width seen by H at the n-th vibrational level. We 
obtain T(^2) = 5.1 x 10"^ and T{Es) = 0.24. H trans- 
fer will occur as a result of climbing a n-step ladder and 
tunneling from this energy level with a combined transfer 
rate RnT{En). Since R2T{E2) > RsTiEs) > R^, at the 
experimental currents of 1-150 pA, n = 2 is the optimum 
level for the combined mechanism. Figure |4] shows R2 
for It = 100 pA. Values at different It are estimated by 
noting that Rn scales as IJ^ in the weak inelastic limit. 
This scaling also shows that H transfer from the sec- 
ond and third vibrational states becomes comparable at 
larger currents {It > 1 nA). Importantly, the similar rates 
obtained with positive and negative V are consistent with 
the experimental insensitivity to the bias polarity. Fur- 
thermore, we observe in Fig. [4] a tendency to yield larger 
rates at F ^ 0.2 V, consistent with the slightly wider 
stripes reported by Sykes et al. at F > 0. This effect is 
a direct consequence of the distinct excitation dynamics 
of holes and electrons. More precisely, if the focusing of 
d carriers had not been considered, the huge asymmetry 
in F|'^ (see inset in Fig. [3| would have favored negative 
voltages over positive ones in Rn {g- g. > 1 order of 
magnitude in R2). 

In conclusion, we unravel the mechanism behind H dif- 
fusion in bulk Pd assisted by inelastic interaction with 
STM-injected hot electrons. Experiments reported for 
this system exhibit comparable H transfer for y > and 

V < 0. This is a counter-intuitive result considering the 
strong difference between electron and hole band struc- 
tures in Pd. We find that both observations are recon- 
ciled by the fact that carriers from Pd d bands do not 
assist H diffusion, as their propagation is restricted to 
narrowly focused cones. The H vibrational excitation is 
only driven by the remaining s carriers, more symmetri- 
cally distributed around Ep- Moreover, we find that at 
the experimental currents of 1-150 pA quantum tunnel- 
ing is crucial to explain H diffusion (without it, the rates 
would be at least 2 orders of magnitude lower). Our 
model allows us to quantify the dependence of the diffu- 
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FIG. 4. (Color online) H transfer rate from the third excited 
vibrational state as a function of the bias voltage for It — 
100 pA. Inset: number of available initial states per Pd atom. 



sion rates on the experimental It and V at low tempera- 
ture. Since this mechanism is based in general concepts 
of solid-state theory i.e., the electronic band structure 
in a periodic potential, the propagation according to re- 
lated Green functions, the e-ph interaction and tunneling 
of atoms at low temperatures, we believe our conclusions 
should apply to other transition metals. 
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